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Abstract 

We perform extensive simulations of a 2D LJ glass subjected to quasi-static shear deformation 
at T = 0. We analyze the distribution of non-affine displacements in terms of contributions of 
plastic, irreversible events, and elastic, reversible motions. From this, we extract information 
about correlations between plastic events and about the elastic non-affine noise. Moreover, we 
find that non-affine motion is essentially diffusive, with a clearly size-dependent diffusion constant. 
These results, supplemented by close inspection of the evolving patterns of the non-affine tangent 
displacement field, lead us to propose a phenomenology of plasticity in such amorphous media. 
It can be schematized in terms of elastic loading and irreversible flips of small, randomly located 
shear transformation zones, elastically coupled via their quadrupolar fields. 
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I. INTRODUCTION 



Plastic deformation of amorphous solids and, more generally, of jammed disordered media 
(foams, confined granular media, colloidal glasses,. . . ) has been intensively studied in the 
past two decades. General agreement is now gradually emerging about the nature of the 
elementary dissipative events in these highly multistable systems. They consist in sudden 
rearrangements of small clusters comprising a few basic structural units, such as Ti events in 
dry foams. In the case of glasses, where they cannot be observed directly, progress has come, 
following the pioneer work of Argon and collaborators, from extensive numerical studies hj, 

nn 

|2, l3|. More recently, simulations performed on model systems - Lennard- Jones (LJ) glasses 
- have proved very helpful to improve our understanding of the effect of topological disorder 
on the elastic as well as plastic shear response of amorphous solids. 

MD simulations are instrumental in elucidating the thermal dependence of the flow stress 
(7(7) in the high strain rate (7 >> 1) regime. However, such conditions (finite T, high 7) 
"blur" the microscopic motion, making it difficult to characterize precisely the elementary 
events, which are the building blocks on which constitutive laws should be based. For this 
purpose, a second class of numerical works have focussed on the athermal (T = 0), quasi- 
static (7 — > 0) regime. In this later regime, hereafter abbreviated as AQS, when a sample is 
sheared at constant rate, the stress-strain cr(7) curve exhibits (see Figure [1]) elastic branches 
interrupted by discontinuous drops Act which are the signature of the dissipative events. 
Beyond an initial transient, o"(7) fluctuates about an average value a, which is identified 
with the yield stress ay = \ima{'y). The distribution of stress drops is broad and system- 
size dependent. In their study on a 2D LJ glass, Maloney and Lemaitre 4] were able to 
analyze them in terms of cascades of elementary events, which we will term "fiips". Each 
such fiip involves both the strong rearrangement of a small cluster ( ~ a few atoms), and 
the appearance of an associated quadrupolar elastic field. This result substantiates the 



6| shear transformations 



representation of elementary processes as Eshelby-like 

However, it appears desirable to analyze AQS simulations in more detail, in order to 
shed light upon debated questions concerning phenomenologies based on the notion of shear 
transformation zones (STZ) s], [s], 1^. Namely: 

- Can the fiipping clusters be associated with regions of the disordered solid which re- 
tain their identity over a finite range of elastic loading before they reach their instability 
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threshold? 

- If so, do the simulations give information about the response of a zone to the elastic 
field generated by the flip of another zone? 

- Can one evaluate the relative importance of the dynamic noise resulting from this 
mechanism as compared with disorder-induced fluctuations of the non-affine elastic field? 

In order to address these issues, we extend in this article a recent study on a 2D LJ 



glass 



10| , by Tanguy et al, which confirms that plastic flow is spatially heterogeneous. They 



claim that one should distinguish between two types of plastic events: strongly localized ones 
occurring during the initial loading phase (a << <ty), and non-localized ones which they 
term "non-permanent shear bands". They also investigate atomic motion in the direction 
transverse to the plastic flow. While they find it to be diffusive at long times, they conclude 
to its hyperdiffusive nature on short time intervals. 

Transverse displacements are purely non-affine and, as indicated by the jagged shape 
of the cr(7) curve, consist of a succession of possibly noisy elastic episodes interspersed 
with sudden atomic rearrangements associated with the plastic events, hence the interest 
of a detailed analysis of their dynamics. With this remark in mind, we revisit in Section 
II the analysis of transverse particle motion, now explicitly separating elastic and plastic 
contributions. We thus extract information about (i) correlations between plastic events (ii) 
the elastic contribution to the non-affine noise invoked in the phenomenological STZ model 
of Falk and Langer 9|]. Moreover, we find it inappropriate to qualify global transverse 
particle motion as hyperdiffusive at short times. Indeed, the effective diffusion coefficient 
D{A'~f) = (A?/^)/A7 smoothly increases with A7 from a finite short time value Do fo the 
asymptotic value Doo- Noticeably, the whole D{A'-f) curve is system size dependent - a result 
which appears consistent with the analysis of plastic events in terms of flip cascades. We 
also show that the contribution of plastic events to transverse diffusion dominates markedly 
over the effect of disorder-induced non-affine elastic fluctuations. 

In Section III, on the basis of close inspection of the evolution of the spatial structure 
of the inflnitesimal non-affine field, we show that our results can be interpreted in terms of 
the elastic loading of zones driven by shear, which soften gradually as they approach their 
spinodal limit. As they approach this threshold, these zones give rise to quadrupolar fields 
of growing amplitude. The zones thus identified can be traced back over shear intervals 
substantially larger than the average interval between plastic events, the effect of which is 



3 



thus easily observable. 

These observations support a description in terms of elastically loaded zones, and 
point towards the relevance of the dynamical noise generated by plastic flips themselves. 

ings should be responsible for the "autocatalytic 



The associated inter-zone elastic coup 
avalanches" j^, [l^ or flip cascades {4, 



13| constituting the system-size dependent plastic 
events, which we believe to be precisely the non-permanent shear bands invoked in ref . . 
We attribute the stronger localization of initial events to the smaller density of nearly un- 
stable zones in the as-quenched or weakly stressed samples. 

This empirical study thus leads to the emergence of a phenomenology of the non-affine 
shear response summarized in Section IV. While supporting the concept of shear transfor- 
mation zones (STZ) or, equivalently, elastically loaded traps (SGR), it diverges from these 
models about two of their basic assumptions, namely independence of elementary events 
and the nature of noise. We think that it should be of use for further developments in the 
modelization of plasticity of amorphous media. 

II. TRANSVERSE PARTICLE MOTION 
A. 2D simulations in the AQS regime 

We use here the same binary LJ mixture as that of ref. ^, namely large (L) and small(S) 
particle radii and numbers are: = 0.5, = 0.3, and Nl = Ns{l + V5)/4. These 
values ensure that no crystallization occurs at low temperature. Simple shear deformation 
is imposed using Lees-Edwards boundary conditions. We study systems of three different 
sizes L X L, with L = 10, 20, 40. 

The quasi-static regime corresponds to the limit where the external time scale 7"^ is much 
larger than that of internal relaxation processes. The system, starting from local equilibrium, 
thus follows adiabatically the shear-induced evolution of the corresponding energy minimum 
up to the spinodal limit where this minimum disappears and the system jumps into another 
local minimum of lower energy. The simulation proceeds as follows: a small increment S of 
homogeneous shear strain is imposed, then energy is minimized using a conjugate gradient 
algorithm with a stringent convergence criterion (see ^ for details). We choose S = 10~^, 
small enough to ensure that, for our system sizes, all elastic branches are well resolved. A 
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typical (7(7) curve is displayed on Figure [H 

Starting from an initial quench, we explore the shear range 7 < 2. In order to characterize 
the stationary state, we only retain data for 7 > 0.1, which ensures that initial transients 
are discarded, thus making details of the quenching protocol immaterial. We have been able 
to collect data on 100, 50 and 20 systems of respective sizes L = 10, 20, 40. 

Following ref.j4|, we make extensive use of the so-called non-affine tangent field {ui}, de- 
fined as the difference between the linear response of particle displacements to an increment 
of homogeneous strain and the corresponding homogeneous field. It is well-defined, and is 
computed, everywhere along each elastic branch. 

B. Statistics of transverse displacements 

In plastically deforming amorphous systems, non-affine displacement fields contain infor- 
mation about both, departures from standard continuum elastic behavior and the nature 
of plastic events. In our simple shear geometry, transverse particle displacements Ai/i are 
purely non-affine, while longitudinal ones mix affine and non-affine contributions. We thus 
focus on the normalized distributions P(Ay, A7) of the Ay^'s for a fixed strain interval 
A7. The statistical ensemble is built by sampling, for all initial configurations, all the 
^Vi = Viilo + A7) - ?/i(7o) at each step 70 (0.1 < 70 < 2 - A7V 

Our results agree qualitatively with those of Tanguy et al lOj] . Namely (see Figure [2]) 
P{Ay, 6) has a quasi-gaussian center and exhibits, beyond Ay ~ 0.2, an exponential tail. Its 
fine structure shows more clearly in the log-lin representation of the distribution -P(C) of the 
scale variable ( = log(A?/). As displayed on Figure [3l for small strain intervals (A7/5 ~ a 
few units), the corresponding curves exhibit a peak at small Ay in the 10~^ range, together 
with a broad hump for larger values. As A7 increases, the peak shifts to the right and the 
hump amplitude grows, though without noticeable horizontal shift, until both merge, for 
A7 ^ 10-2. 

Can we interpret this structure and its evolution in the light of the succession of elastic 
and plastic episodes which reflect into the sawtooth shape of the cr(7) response? Clearly, 
particle motion consists of a series of continuous trajectories interrupted by sudden jumps 
associated, respectively, with elastic episodes and plastic events. Since 6 is our numerical 
strain resolution, each interval of length 6 is, within our accuracy, defined as either purely 
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elastic or purely plastic. So, for A7 = S, we can decompose P (and likewise P) as 



P(Ay, 5) = api{5)Upi{Ay, 5) + (1 - api{5))U,i{Ay, 5) (1) 

where ftp; (5) is the fraction of "plastic intervals" in our ensemble, and Upi (resp. He/) are 
the normalized distributions associated with the plastic (resp. elastic) sub-ensemble. It now 
clearly appears (see Figure HI) that the small Ay behavior of P{Ay,6) results entirely from 
elastic motion, while its hump and large-Ay tail are due to plastic events. 

On this basis, we are now able to understand the evolution of P with A7. We extend 
the above decomposition to A7 > 6, now defining as "plastic" any interval containing at 
least one plastic event. As A7 increases while remaining small with respect to the average 
length {6api{6)~^)) of an elastic branch, most plastic intervals in general still contain only 
one plastic event, so that their fraction api{Aj) ^ {Aj/6)api{6), while Upi remains quasi 
unchanged, because the scale of plastic slips is large compared with that of elastic ones. 
This explains the evolution (Figure H]) of the plastic hump which, as A7 increases, consists 
primarily in an upward shift with little change in shape. 

In this regime, which corresponds to A7/5 << api{6)^^ ( ~ 100 for the 20 x 20 system), 
the fraction aei = (1 — api) of purely elastic intervals remains nearly constant (e.g., for 
L = 20 and A7 = 105, 6aei/aei ~ 0.1), so that the variations with A7 of the elastic part of 
P, Pel = aei{Aj)Ilei, directly reveal those of Ilez. Inspection of the numerical data suggests 
(see also Figure H]) that the rightward shifts of Pei{C) ^-re roughly equal to log(A7). The plot 
(Figure [5]) of Pei vs \og{Ay/Aj)) shows that the collapse is indeed excellent in the low and 
middle ranges of {Ay/A'-f) values, but fails in the large Ay tails. 

This collapse means that during the purely elastic episodes separating plastic events non- 
affine displacements are essentially convective: the 7-derivative of the yi field is nearly 



constant over most of an elastic episode 
field juj j, although spatially disordered 
ref. 



his amounts to stating that the tangent non-affine 



llj, remains quasi-quenched. Now, we know from 
141 ] that this cannot be true when the system comes close below a spinodal threshold 7c. 
Indeed, in these near-critical regions, a gradually softening elastic mode develops, leading to 
a (70 — 7)"^'^^ divergence of {ui}. Clearly, it is the contribution of these near-critical softened 
configurations which explains the large- Ay tails of fig/. The tail deflation (Figure [5]) results 
from the fact that, as the interval A7 increases, the weight of these soft configurations is 
gradually transferred from Yl^i to 11^;. 
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This analysis clarifies the physical significance of the shape of P{Ay, A7) in the moder- 
ate A7 range (where its non-Gaussian parameter remains large 1^): (i) Its quasi-gaussian 



center 



111 ] results from small-scale non-affine displacements accumulated along purely elas- 



tic segments, during which particle trajectories are essentially convected. (ii) Its quasi- 
exponential tail arises from plastic jumps during plastic events. 

For larger ^ api{6)~^, however, the decomposition of equation ([1]) loses physical 

content since the statistical weight of elastic intervals vanishes while Ilj,; mixes both elastic 
and plastic contributions. In order to circumvent this limitation, we concentrate in the 
following on the evolution with A7 of the second moment (Ay^). 



C. Transverse particle diffusion 

In order to elucidate the nature of the transverse particule dynamics, we have computed 
from our data the space and ensemble average (Ay^) for increasing values of A7. On 
Figure [6] we plot the effective diffusion coefficient D{A'~f) = {Ay'^)/A'-f. For our three 
system sizes, D{A'-f) exhibits the same qualitative features, namely it increases from a finite 
value Dq = lim D, and saturates at a finite value D^o for large A7 0.5). However, D is 
strikingly system-size dependent, D^^\A'j) increasing with L for all A7's. 

Again, in order to shed light on the origin of this behavior, let us separate explicitly the 
elastic and plastic contributions to non-affine particle motion. Indeed, we can write: 

^=U^y{l)+Y.y^''Sil-la) (2) 

where a labels the strain values where irreversible events (avalanches) occur, so that : 

Z)(A7) = Z)ee(A7) + I)ep(A7) + DppiA^) (3) 

with 



Dee{A-f) = 




and -Dgp, Dpp the corresponding cross-correlated and plastic-plastic contributions. These 
three functions carry the information about self and cross correlations of elastic and plastic 
non-affine displacements. 

Let us first consider the A7 — > limit. A straightforward asymptotic analysis taking into 
account the square-root divergence of the tangent field {ui(7)} at spinodal points leads to: 

Dep(A7) = 0(A7^/^) (5) 
7 



Dee(A7) = 0(-A7Log(A7)) (6) 

So, their contribution to D(A7) vanishes, and the finite value of Dq originates only from 
avalanches: 

Do=]hn^D,, = f {{Ytf) (7) 



with 

/ = hm 



the average avalanche frequency, and {{Y^^y)) the variance of transverse displacements in a 
single plastic event. 

It was shown in 0, [l^ that / increases with size roughly as the lateral system size L. 
Using the measured frequency, we find that the variance = -Dq/"^ decreases slowly 

with size: for L = 10,20,40, we get lO^iiYi^f)) = 37,26,20. We will come back to this 
point later. 

Figure [7] shows the decomposition (eq.Q) for L = 20. It is seen that, for all A7, 
Dpp remains the dominant contribution. (A7) grows, then saturates. Its growth range 
Tpp ~ 0.25 shows no clear size dependence. This behavior entails that plastic events are 
correlated over a finite A7 range, which we understand as measuring a typical shear range 
over which the structural features involved in plastic events retain their identity. 

For all sizes, Dep{A'-f), though much smaller than Dpp, exhibits an analogous behavior 
with a similar correlation range. This indicates that the above-mentioned persistent struc- 

n 

tures dominate the non-affine elastic response. Moreover the near square-root behavior pj[ 
of Dep for small A7 signals that the field {ui} in a given near-critical region is strongly 
correlated with the displacements associated with the subsequent avalanche. 

Finally, Z^ee also rises and saturates, but more slowly than D^p and Dpp. From the data 
for the three system sizes, we evaluate its correlation range to be Fee ~ 1. That is, a single 
renewal, after A7 ~ Tpp, of the "active structures" is not sufficient for the non-affine field 
to fully decorrelate. This we take as a hint of the fact that active structures occupy on 
average a fraction only of the total system "volume". With this interpretation, we evaluate 
this fraction as Fpp/Fee ~ 1/4. 
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III. ZONE EMERGENCE, FLIPS AND ELASTIC COUPLINGS 



Two phenomenological models of plasticity of jammed disordered media - namely the 
STZ theory of Falk and Langer j^, Q, and the soft glass rheology (SGR) of SoUich et 
al Q - have been proposed recently. Both describe, explicitly (STZ) or implicitly 
(SGR), plastic events as transitions concerning small regions, modelized as either Eshelby- 
like transformations or jumps out of traps in the energy landscape. These flips are viewed 
as independent, hence individual, and governed by the combined effect of external drive and 
of a thermal-like noise. 

A question then immediately arises : are our above results consistent with the basic 
ingredients of these models? We now try to answer it by careful inspection of the evolution 
of the spatial structure of the non-affine tangent field {ui}. 



A. Identification of active structures 



Let us first recall that {ui} can be written as 



Uia = (9) 

- ~ (10) 

where label particles and (a,/5) cartesian components. U is the total energy of the LJ 
system, and H the associated hessian matrix. 3^7 is the field of particle force increments 
generated by an infinitesimal increment of homogeneous shear, {ui} can be decomposed on 
the eigenmodes ^'n of H (the phonon modes). Upon approaching a spinodal point, a single 
one, ^'i, softens critically Q, Q. So, in near critical regions, {ui} exhibits a square-root 
divergence and is dominated by its projection ui on ^'i. This property enables us to clearly 
characterize the active structures close to the onset of plastic events. 

The features we identify as characteristic of the evolution of {u} are exemplified on 
Figures 8 and 9. In these figures, for each value of 7 we decompose the total non-affine field 
as: 

U = Ui U2 u (11) 

with ui^2 its projections on the two modes ^'1^2 with lowest eigenvalues (Ai < A2). 
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We observe that, in all cases, the active structures emerge out of a quasi-random small 
background as one or several localized zones characterized by strong quadrupolar contribu- 
tions to u. These quadrupoles are approximately aligned with the principal directions of 
the homogeneous strain. As zones approach instability, they soften, and the corresponding 
growing part of u concentrates into the lowest mode component Ui. 

After identifying them in prespinodal regions, we are able to trace them back over sizeable 
7-ranges before their amplitude has decreased enough for them to gradually merge into the 
global disordered background structures. In some cases, this range of visibility extends across 
one or more plastic events. A typical sequence illustrating this behavior, shown on Figure [HI 
extends over a range of A7 ~ 3%, to be compared with the average length A7 ~ 1% of 
elastic episodes for this system size. 



B. Zone flips 

Each discontinuous drop on the cr{'-f) curve is associated with the sudden disappearance 
of one or more of these soft zones (Figures [8] and [9]). Noticeably, the resulting changes in 
{ui} remain quite localized, leaving most of the other prominent structures of the non-affine 
pattern essentially unchanged. In particular, zones which were already clearly visible are 
commonly seen to survive the plastic event (see Figure [H]). 

We have shown in Section II. 2 that particle jumps {ARf} in plastic events are well 
correlated with the precritical {ui} structure. This entails that the flip of a zone Z is 



associated primarily with a quadrupolar displacement field of finite anip 

3, 



itude centered on Z. 
8| and explicited by 



This supports the representation, proposed by Argon and coworkers 

Picard et al 3] , of the elementary dissipative process as the Eshelby-like shear transformation 
of a self generated inclusion involving a few particles, which is also the basis of the STZ 
theory. 

This representation is also consistent with the main features of the distribution Ilpi[Ay, S) 
of particle displacements induced by plastic events (see § II. B). Indeed, the plot on Figure [TOl 
shows that that Upi becomes exponential for Ay ^ 0.2. In the Eshelby picture, we expect 
large Ay to correspond to displacements within the transforming zone(s), while smaller Ay 
are associated with particles sitting in the surrounding elastic medium. In this picture, 
we interpret the exponential tail of lie;, hence of the full distribution P, as reflecting the 
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diversity of intrazone structures. We thus expect that, upon varying the system size, the 
logarithmic slope of the tail should remain constant. As shown on Figure [101 this prediction 
is very nicely verified for our three L values. Moreover, it appears that the tail amplitude 
decreases with increasing L. In our interpretation, the statistical weight of the exponential 
tail is controlled by the volume fraction of zone cores involved in plastic events. Since the 
average size of the avalanches constituting plastic events scales roughly as L, we expect the 
corresponding core volume fraction to decrease with size, in agreement with the observed 
behavior. 

So, we interpret the existence of an exponential tail in P{Ay) as a consequence of struc- 
tural disorder within transforming zones, and by no means as a signature of avalanche 
behavior, which only affects the tail amplitude. 



C. Elastic couplings 

Consider the sequence shown on Figure [9l For 7 = 0.1150 (Figure [9]-a) two zones are 
clearly discernible in the u pattern. At 7 = 0.1160 (Figure [9]-b) their amplitude has grown, 
and they appear in the projection U2 of u on the next-to-lowest mode. At 7 = 0.1161 
(Figure [9]-c) they are slightly softer and have invaded the projection Ui on the lowest mode. 
Note that, in these last two cases, they appear in the non-trivial soft mode as connected by 
"flow lines" which reproduce the most prominent vortex structures first described by Tanguy 
et al 20|. This is particularly clear when two relatively distant zones soften simultaneously 
as shown on Figure [TTl 

From these and many similar observations (see Figure [8]) we deduce that : 

• Soft zones are coupled elastically via their quadrupolar fields. 

• It is the associated flow lines which form the conspicuous vortex-like patterns character- 
istic of non-affine fields in amorphous solids. That such patterns are akin to incompressible 
flows results from the fact that these systems are much more compliant in shear than in 
compression. For our 2D LJ glass ^: /i = 39, K = 236. 

Quadrupolar couplings result, in the homogeneous elastic continuum approximation 
for the background medium, from a stress field ~ cos(46')/i?^ (with R the length and 6 
the orientation with respect to x of the interzone vector). Notice that, in Figure in ui 
the fields of the two zones Zi, Z2, which are quasi vertically aligned, combine positively. 
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We have checked that, in consistence with the above remark, zone pairs lying in sectors 
corresponding to negative couphngs combine negatively in the phonons ^'i and ^'2: in 
this case, one of the two zones appears with an opposite sign, i.e. with reversed flow lines. 
These observations conflrm the relevance of inter-zone elastic couplings. 

It is intuitively clear that conflgurations with strongly coupled soft zones are good can- 
didates for simultaneous flips, as illustrated on Figure [8] and |9l We observe a number of 
multi-zone flips, involving the disappearance of at least one diverging quadrupole, as well 
as that of some less visible zones. Moreover, plastic events alter the amplitude of the near- 
critical surviving zones in u (see Figure [S]) in a way which appears roughly consistent with 
quadrupolar couplings. 

IV. A PHENOMENOLOGY OF SHEAR ZONES 

The various pieces of information which we have gathered and presented above can now 
be organized into a rather detailed phenomenology of plasticity in our system in the AQS 
regime, which can be schematized as follows. 

(i) Structural disorder gives rise, in such a glassy system, to the existence of strong inho- 
mogeneities, or zones. These zones can be viewed as small inclusions a la Eshelby, embedded 
in a quasi-homogeneous elastic background, and plastically transformable under shear. They 
are advected by external shear toward their instability threshold. Upon approaching it, a 
zone softens, then flips at its spinodal and disappears. 

The correlation range Tpp between plastic events (see Section II. 2) provides an evaluation 
of the amount of strain necessary to fully renew the population of zones. Here Tpp ~ 0.25. 

One renewal is not suflicient, however to decorrelate fully the non affine fleld. This 
suggests that zones occupy a fraction only of the system, of order Tpp/Tee, here ~ 1/4. 

(ii) Zone elastic softening and plastic flipping are both associated with quadrupolar com- 
ponents in the non-affine displacement fleld which give rise to inter-zone elastic couplings. 
Due to the long range of elastic flelds, a zone flip thus alters the local strain level at any 
other zone site in the system, the resulting 7-shift depending, in amplitude and sign, on the 
relative position of source and target. These signals have two types of effects. 
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— A flip signal may shift the strain level of some zones beyond their threshold 7c, thus 
triggering their flip and initiating a cascade. In the QS regime, where acoustic delays are 
neglected, such cascades are instantaneous. 

— For the non-flipping zones, the elastic signals resulting from flips constitute a dynamical 
noise, acting in parallel with the externally imposed advection, whose frequency scale is 
proportional to the strain rate 7. 

So, in agreement with the previous proposition of Bulatov and Argon [s], Q], also 
underlying recent models by Picard et al 2l| and by Baret et al 22|], we conclude that 
elastic couplings between zones play an essential role, by inducing avalanches and generating 
the dominant contribution to the disorder-induced noise. At T = and moderate strain 
rates, it is this dynamical elastic noise which, in addition to strain advection, must appear 
in models of plasticity of jammed media. In a forthcoming article, we will propose such 
a schematic model, and show that it accounts for a system-size scaling behavior of avalanches. 



Within the foregoing picture, we interpret the two kind of events invoked by Tanguy et 
al [l^ as plastic events of different sizes. Namely, their transient shear bands exhibit the 
characteristics expected for large avalanches of the type described above, their directionality 
resulting from that of quadrupolar couplings in their system with two rigid walls. Their 
local events correspond to the small stress drops appearing during initial transients (see 
Figure [1]): at low stress levels, before advection has been able to significantly feed the 
near-threshold population, we indeed understand that avalanche triggering is unlikely, 
hence that single flip events are the rule. 



Let us finally stress that the above zone phenomenology remains very schematic. Indeed, 
exhaustive inspection of the evolution of u reveals that zone life is somewhat more eventful 
than our simplified description suggests. Rather frequently, we see a given zone undergoing 
a few flips before it disappears. Such zones can be termed multi-state. In other instances, a 
barely visible zone emerges fast enough to "overtake" previously more visible, softer, ones. 
That is, zone moduli and, very likely, thresholds are not unique, but distributed about an 
average. These remarks point toward the interest of pursuing extensive characterization of 
elastic heterogeneity in jammed systems, in particular via studies of coarse grained elastic 
moduli. 
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On a more speculative level, we would like to raise an important question: are the 
irreversible transformations identified here under shear related to the dynamic hetero- 



geneities [23 



24 



25 



26(1 characteristic of glassy dynamics near and below Tgl Indeed, 



we share the opinion, formulated long ago by Goldstein [^^j, that at finite temperature "re- 
arrangements are of course occurring all the time in the absence of an external stress; the 
external stress, by biasing them, reveals their existence" . 

Our active zones are primarily sensitive to shear. This, we think, must be put together 
with recent results by Widmer-Cooper and Harrowell 28|, |29l]. These authors show that 
dynamic propensity in a 2D LJ glass is uncorrelated with free volume, which we understand 
to mean that their local rearranging structures are only wekly coupled to compression. This 
leads us to suggest that their observation that zones of high propensity have large " Debye- 
Waller factors" amounts to identifying them as soft zones sensitive to shear - a speculation 
which will demand extensive future investigation. 
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FIG. 1: Stress-strain curve for a 20 x 20 system. 
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FIG. 2: The distribution P{Ay, 5) of transverse particle displacements for the elementary strain 
step 6 and for system size L = 20. Insert: Blow-up of the small Ay region. 
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FIG. 3: The distribution -P(C; A7) of the scale variable C = log(Ay) for increasing values of 
n = A7/(5 = 1,4, 16, 64, 256. System size: L = 20. 




FIG. 4: Decomposition of the total distributions A7) (thick solid lines), for n = Aj/6 = 
1,8,64,512. For increasing n's, the maximum of the distribution shifts rightwards. Thin solid 
lines: contribution of plastic events; thin dashed lines: contribution of elastic branches (see text). 
System size L = 20. 
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FIG. 5: The function plotted versus Ay/A^ for n = A^/5 = 1,2,4,8,16,32,64. The arrow 
indicates the direction of increasing n's. System size L = 20. 
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FIG. 6: From bottom to top: the instantaneous diffusion constant, D{A'y) = {Ay^) / A'y versus 
A7, for increasing system sizes L = 10,20,40 (and respectively 100, 50, and 20 samples). Error 
bars result from a standard Student analysis of the data. 
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FIG. 7: The instantaneous diffusion constant, D{A'y) = (A7/^)/A7 for system size L = 20 (thick 
sohd hne) and its decomposition into D^e (dotted hne), D^p (dashed hne) and Dpp (thin sohd hne). 
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FIG. 8: System size: L = 20. Each line corresponds to a single strain value 7. The four frames 
present, from left to right, the non-affine tangent field u and its three components ui, U2, u as 
defined by equation ([TT]) . 

(a) : 7 = 0.0556: a quadrupole, clearly visible in the upper right corner of both u and the soft 
mode component ui, signals a near-critical zone Z, which flips at 7c G [0.0557,0.0558]. 

(b) : 7 = 0.0558: Z has just flipped and disappeared, pushing another zone Z' (lowest right 
corners of u and ui) closer to its threshold. Note that Z' was already discernable (line a) in both 
u and u before the event. ^'^ 

(c-d): Zone Z, now indicated by circles, can be traced hack to 7 = 0.0269 (line c, frame u). It 






FIG. 9: From a to d: 7 = 0.1150, 0.1160, 0.1161, and 0.1166. Two vertically aligned zones flip 
together in the interval 0.01165 — 66 and disappear. As strain increases from 7 = 0.1150, they soften 
and sequentially appear in U2, then in ui after a clear level crossing in the interval 0.01160 — 61. 
The connecting flow lines are a signature of their elastic coupling. 
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FIG. 10: The distribution Ilpi{Ay, 5) of particle displacements induced by plastic events for L = 
10,20,40. 




FIG. 11: Two distant zones soften simultaneously and appear in ui. The characteristic vortex-like 
structure signals their elastic coupling. 
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